Resummation of Nonalternating Divergent Perturbative Expansions 
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A method for the resummation of nonalternating divergent perturbation series is described. The 
procedure constitutes a generalization of the Borel-Pade method. Of crucial importance is a special 
integration contour in the complex plane. Nonperturbative imaginary contributions can be inferred 
from the purely real perturbative coefficients. A connection is drawn from the quantum field theoretic 
problem of resummation to divergent perturbative expansions in other areas of physics. 
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In view of the probable divergence of quantum field 
theory in higher order 0J^] , the resummation of the per- 
turbation series is necessary for obtaining finite answers 
to physical problems. While the divergent expansions 
probably constitute asymptotic series H, it is unclear 
whether unique answers can be inferred from perturba- 
tion theory [QJ^]. Significant problems in the resumma- 
tion are caused by infrared (IR) renormalons. These are 
contributions corresponding to nonalternating divergent 
perturbation scries. The IR renormalons are responsi- 
ble for the Borel-nonsummability of a number of field 
theories including quantum chromodynamics (QCD) and 
quantum electrodynamics (QED) [|4|||. 

Here I advocate a modification of the resummation 
method proposed in for nonalternating divergent 

perturbation series. The method starts with a given in- 
put series, 



/(<?) 



E 



c„ > 0. 
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where g is the coupling parameter and the perturbative 
coefficients c n are expected to diverge as follows ||, 



K 



S" 
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with K, 7 and S being constant. The Borel transform 
/b of the perturbation series (111) 



(3) 



has a finite radius of convergence about the origin. For 
the evaluation of the Borel integral, /b(<7) has to be 
continued analytically beyond the radius of convergence. 
Strictly speaking, this analytic continuation has to be 
done on the branch cut in view of the nonalternating 
character of the series (Q). This requirement can be re- 
laxed slightly by performing the analytic continuation 
into regions where g acquires at least an infinitesimal 
imaginary part g — ► g ± i e. In this case, the analytic 
continuation can be achieved by evaluating Pade approx- 
imants M . The first n+l terms of the Borel transformed 



series (||) can be used to construct a diagonal or off- 
diagonal Pade approximant (for the notation see ) 



r n (z) = 



[n/2]/|[(n + l)/2] 



(4) 
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where JxJ denotes the integral part of x. The resum- 
mation is accomplished by constructing the sequence of 
transforms {Tf n (g)}^— where 



Tf n (g) 



dt exp(-t)V n (gt) , 



(5) 



and the integration contour Cj (j = — 1,0, +1) is as 
shown in Fig. [j] (for j = — 1 and j = +1). The result 
obtained along C_i is the complex conjugate of the re- 
sult along C+i. The arithmetic mean of the results of 
the integrations along C-\ and C+i is associated with 
Cq- Therefore, the result along Cq is real rather than 
complex. The limit of the sequence {T f n (g)}^Lo (pro- 
vided it exists), 



lim Tf n (g) = f(g), 



(6) 



is a plausible complete nonperturbative result inferred 
from the perturbative expansion ([l]) . Which of the con- 
tours Cj (j = —1,0, +1) is chosen, has to be decided on 
the basis of additional considerations which do not follow 
from perturbation theory alone. 

The zeros of the denominator polynomial of the Pade 
approximant [see Eq. (|4|)] correspond to the poles of the 
integrand in Eq. (|^). Denote by t the integration vari- 
able for the evaluation of the generalized Borel integral 
in Eq. (||), then the poles lie at t = Zi (where the in- 
dex i numbers the poles) along the positive real axis 
(Imzi = 0) and in the complex plane (Imz; ^ 0). The 
poles lying on the positive real axis are treated as half- 
poles encircled in the mathematically positive sense for 
C_i and as half-poles encircled in the mathematically 
negative sense for C+i. The contour C_i encircles all 
poles at t = Zi in the lower right quadrant of the com- 
plex plane (Rez^ > 0, Imzj < 0) in the positive sense 
(see Fig. [j]). The contribution of these poles should be 
added to the final result. The contour C+i is understood 
to encircle all poles in the upper right quadrant of the 
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complex plane in the mathematically negative sense. In 
general, the integrations along C_i and C+i lead to a 
nonvanishing imaginary part in the final result for f(g) 
[see Eq. @], although all the perturbative coefficients 
c n are by assumption real and positive [see Eq. O)]. It 
might be interesting to note that, as with any complex 
integration, it is permissible to deform the integration 
contours shown in Fig. [l] in accord with the Cauchy The- 
orem as long as all pole contributions are properly taken 
into account. 

FIG. 1. Integration contours for the evaluation of the 
generalized Borel integral in Eq. (0). 
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This paper represents a continuation of previous work 
on the subject ^[^J]- The resummation method defined 
in Eqs. (0)-(||) differs from || in the combination of Borel 
and Pade techniques and, if compared to the remarkable 
investigations in [Q,^) on the resummation of QCD per- 
turbation series, in the integration contour used for the 
evaluation of the generalized Borel integral. It is argued 
here that, when the Borel transform (|3|) is analytically 
continued with Pade approximants (^), the contribution 
of poles lying off the positive real axis has to be taken into 
account in order to obtain consistent results in the resum- 
mation (see Fig. p. In j7|,[)| it is argued that the Borel 
integral should be evaluated by principal value. It could 
appear that the Co contour corresponds to the principal- 
value prescription. However, this is not necessarily the 
case, if there are poles present which lie off the positive 
real axis (i.e., at t = Zi with Rezj > 0, Ivazi ^ 0). The 
contribution of these poles not only modifies the imagi- 
nary, but also the real part of the final nonperturbative 
result. Of course, when there are no poles lying off the 
positive real axis, as it is the case for the problems dis- 
cussed in [^U^[ , then the principal- value prescription used 
in j7||| is equivalent to the Co contour. Because the re- 
sult obtained along Co is real, this contour should be used 
whenever the existence of an imaginary part is discour- 



aged by physical reasons. 

It is important to mention that the method presented 
here is not the only prescription currently available for 
the resummation of divergent perturbative expansions in 
quantum field theory. For example, the 8 transforma- 
tion (see Eq. (4) in [pT[ ) is a very useful method for 
the resummation of divergent perturbation series. The 8 
transformation has a number of appealing mathematical 
properties, including rapid and numerically stable con- 
vergence, and it has been shown to yield consistent re- 
sults in many cases, including applications from quantum 
field theory jlTJ and from other areas of physics [|l2| . Be- 
cause the <5 transformation fulfills an accuracy-through- 
order relation (see Eq. (9) in jnj), it can be used to pre- 
dict perturbative coefficients. The 8 transformation is 
primarily useful for alternating series. It fails, in general, 
in the resummation of the nonalternating series discussed 
here. The 8 transformation and the resummation method 
introduced here complement each other. 

Three applications of the resummation method defined 
in Eqs. are considered below: (i) the QED ef- 

fective action in the presence of a constant background 
electric field, (ii) a mathematical model series which sim- 
ulates the expected large-order behavior of perturbative 
coefficients in quantum field theory, (iii) the perturba- 
tion series for the energy shift of an atomic level in a 
constant background electric field (including the auto- 
ionization width). The nonperturbative imaginary con- 
tributions obtained along C_i and C+i find a natural 
physical interpretation in all cases considered. 

The QED effective action, or vacuum-to-vacuum am- 
plitude, in the presence of a constant background electric 
field has been treated nonperturbatively in and 
the result is proportional to the integral 
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where g^ is a coupling parameter proportional to the 
square of the electric field strength, g& — e 2 E 2 /m*. 
Here, m c is the electron mass, and e is the elemen- 
tary charge. The natural unit system (h = c = 1) is 
used. The imaginary part of <S(<7e) is proportional to the 
electron-positron pair production amplitude per space- 
time interval [there is of course also a muon-antimuon 
pair-production amplitude, obtained by the imaginary 
part of (^) under the replacement m c — > rn M , which is 
not discussed here]. S(<7e) has the following asymptotic 
expansion in the coupling parameter, 



%e) ~ 16 
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4 n |B- 



2n+4| 



(2n + 4) (2n + 3) (2n + 2) 

71 — 7 



<?E +1 
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where #271+4 is a Bernoulli number. In view of the asymp- 
totics 



A n \B 2n+i \ T(2n + 2) 



_2n+4 



1 + 



, n — ► 00 , (9) 
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the perturbative coefficients, which are nonalternating 
in sign, diverge factorially in absolute magnitude. The 
asymptotic series (||) for S(<7e) is taken as the input series 
for the resummation process [Eq. (EJ)], and a sequence 
of transforms TSVi(ffE) is evaluated using the prescrip- 
tion (||). The results have to be compared to the exact 
nonperturbative expression (Q). This is done in Table || 
for gE = 0.05. The partial sums of the asymptotic se- 
ries (B) are listed in the second column. 



mg) ~ E 



has been used as a paradigmatic example for nonalter- 
nating divergent series in the literature P Jl5| , ^6| . This 
series can be resummed by the method (pJ). Moreover, 
this resummation is even exact for all transformation or- 
ders n > 2. This can be seen as follows. The Borel 
transform 



TABLE I. Resummation of the asymptotic series for 
the QED effective action (^) in a constant background 
electric field for <?e = 0.05. Results in the third column 
are obtained by the method indicated in Eq. (jsj) along the 
integration contour C_i. The partial sums in the second 
column are obtained from the asymptotic series (0). 
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2 


0.001 


146 032 


0.001 


144 


848 - 


i7.70 


X 


10- 


17 


3 


0.001 


146 705 


0.001 


146 


639 - 


i8.22 


X 


Hr 


11 


4 


0.001 


146 951 


0.001 


147 


113 - 


i3.54 


X 


nr 


8 


5 


0.001 


147 087 


0.001 


147 


264 - 


i 1.93 


X 


10" 


8 


6 


0.001 


147 195 


0.001 


147 


173 - 


13.15 


X 


10" 


7 


7 


0.001 


147 310 


0.001 


147 


113 - 


i2.58 


X 


nr 


7 


8 


0.001 


147 469 


0.001 


147 


162 - 


i2.30 


X 


10" 


7 


9 


0.001 


147 743 


0.001 


147 


165- 


i2.63 


X 


HT 


7 


10 


0.001 


148 327 


0.001 


147 


144- 


i2.53 


X 


1GT 


7 


11 


0.001 


149 825 


0.001 


147 


157- 


i2.46 


X 


nr 


7 


12 


0.001 


154 375 


0.001 


147 


155 - 


i2.56 


X 


nr 


7 


13 


0.001 


170 560 


0.001 


147 


151 - 


i 2.51 


X 


nr 


7 


14 


0.001 


237 137 


0.001 


147 


156 - 


i 2.51 


X 


10" 


7 


15 


0.001 


550 809 


0.001 


147 


153 - 


i2.53 


X 


1GT 


7 


16 


0.003 


228 880 


0.001 


147 


154 - 


i 2.51 


X 


1GT 


7 


17 


0.013 


345 316 


0.001 


147 


154 - 


i 2.52 


X 


nr 


7 


18 


0.081 


610 937 


0.001 


147 


153 - 


i2.52 


X 


nr 


7 


19 


0.594 


142 371 


0.001 


147 


154 - 


i2.52 


X 


nr 


7 


20 


4.852 


426 276 


0.001 


147 


154 - 


i2.52 


X 


10" 


7 


exact 


0.001 


147 154 


0.001 


147 154 - 


i 2.52 x 10" 
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Numerical results from perturbation theory are nor- 
mally obtained by (optimal) truncation of the perturba- 
tion series. For the example considered, (i) the partial 
sums do not account for the imaginary part and (ii) due 
to the divergence of the perturbative expansion, no im- 
provement in the final result could be obtained by adding 
more than the first 7 perturbative terms, ft requires a 
valid resummation procedure to go beyond the accuracy 
obtainable by optimal truncation of the perturbation se- 
ries. The transforms TS n {g^,) displayed in the third col- 
umn of Table [| apparently converge to the full nonpertur- 
bative result given in Eq. (Q), and the nonperturbative 
imaginary part, which corresponds to the pair production 
amplitude, is reproduced although the input series (|^) 
has purely real perturbative coefficients. 

Two specific mathematical model series are considered 
next. The series 



n=0 



1/(1-5) 



is a geometric series. The summation of geometric se- 
ries inside and outside of the circle of convergence by 
Pade approximants is exact in all transformation orders 
n > 2. So, for all n > 2 the transforms TM n {g) ful- 
fill the equality TN n {g) = -1/g exp(-l/<?) r(0, -1/g) = 
J\f(g), where T(0,x) is the incomplete Gamma function 
(see fL7|| ), and the choice of the contour (C_i or C+i) de- 
termines on which side of the branch cut the incomplete 
Gamma function is evaluated. 
The asymptotic series, 



M(g) 



oo „, 

n=0 



7) 



r(n) 



■n\g 



(10) 



constitutes a more interesting application of the resum- 
mation method than Af(g). On account of the asymp- 
totics, 



T(n + 7 ) 
r(n) 



1 



oil 

n 



oo , 



(11) 



the series M. (g) serves as a model for the expected large- 
order behavior of perturbative coefficients in quantum 
field theory [see Eq. The analytic summation of 

© leads to 



M(g) = T( 1 ) [g 



d_ 



2-^0(1,759) 



(12) 



where the hypergeometric 2-F0 function has a branch cut 
along the positive real axis (see [0). The imaginary 
part of (|1J) for g > as a function of g and 7 is 
im.M(g) = 7T (1 — 37)3 _7_1 exp(— 1/g), where the in- 
tegration is assumed to have been performed along the 
contour C+i. For C_i, the sign of the imaginary part 
is reversed. The numerical example considered here is 
7 = 2.3, g = 0.1. In the Table |l[ numerical results 
are dis play ed for the nth partial sums of the asymptotic 
series ( |l0|) and the transforms TM. n (g) calculated ac- 
cording to Eq. (||) in the range n = 2, . . . , 12. While the 
partial sums eventually diverge, the transforms TM n (g) 
exhibit apparent convergence to about 6 significant fig- 
ures in (n = 12)th transformation order, and the trans- 
forms reproduce the imaginary part although the coeffi- 
cients of the series (uw are all real rather than complex. 
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The integration is performed along the contour C+\. The 
exact result in the last row of Table [n] is obtained from 
Eq. fll^). For the evaluation of the transforms TAi n (g) 
it is crucial to use the contour C+i rather than a con- 
tour infinitesimally above the real axis. For example, in 
order to obtain consistent numerical results, it is neces- 
sary to take into account the poles at t = 9.99 ±i 0.578 in 
(n = ll)th transformation order, encountered in the eval- 
uation of the transform TMn(g) according to Eq. (||), 
and the pole at t = 9.99 ± i 0.495 in (n = 12)th order for 
the evaluation of T M. 12(g)- These poles approximately 
correspond to the triple pole at t — 1/(0.1) = 10 which 
would be expected in the case 7 = 2. 

TABLE II. Resummation of the model series (|l^) for 
7 = 2.3, g = 0.1 by the method indicated in Eq. (^) 
along the integration contour C+i. The partial sums are 
obtained from the asymptotic series jlo|). 



n 


partial sum 


TM n (g) 


2 


0.445 451 


0.393 554 + i 0.373 912 


3 


0.559 685 


0.840 561 + i 0.446 830 


4 


0.640 410 


0.764 942 + i 0.274 640 


5 


0.703 981 


0.765 339 + i 0.218 156 


6 


0.759 669 


0.763 012 + i0.219 638 


7 


0.813 594 


0.762 186 + i0.219 197 


8 


0.870 909 


0.762 196 + i0.219 126 


9 


0.937 322 


0.762 224 + i 0.219 123 


10 


1.020 707 


0.762 225 + i0.219 127 


11 


1.133 528 


0.762 223 +i 0.219 127 


12 


1.297 220 


0.762 223 + i 0.219 127 


exact 


0.762 223 


0.762 223 + i 0.219 127 



When an atom is brought into an electric field, the lev- 
els become unstable against auto-ionization, i.e. the en- 
ergy levels £ acquire a width T (£ — > Rc£ — IF/ 2 where 
r is the width). Perturbation theory cannot account for 
the width. The coefficients are real, not complex (l^] . An 
established method for the determination of the width 
is by numerical diagonalization of the Hamiltonian ma- 
trix |lj] It is argued here that the full complex en- 
ergy eigenvalue, including the width, can also be inferred 
from the divergent perturbation series by the resumma- 
tion method defined in Eqs. (|lj)-(||), where the appropri- 
ate integration contour is C+i. Perturbative coefficients 
for the energy shift in arbitrarily high order can be in- 
ferred from the Eqs. (9,13-15,28-33,59-67,73) in @. 

The symmetry of the problem suggests the introduc- 
tion of the parabolic quantum numbers m, TI2 and m ]2^j 
(the principal quantum number is n = n\ + TI2 + m + 1). 
Here, calculations are performed for the ground state 
with parabolic quantum numbers n\ = 0, ri2 — 0, m = 
and two L shell states, both of which are coherent su- 
perpositions of the 2S and 2P states. One of the L shell 
states investigated here has the parabolic quantum num- 



bers ri\ = 1, 712 = 0, m = 0, and the other L shell 
state has the quantum numbers rii = 0, n<z = 1, m = 0. 
The Stark effect is interesting because, depending on 
the atomic state, the perturbation series are either com- 
pletely nonalternating in sign (e.g., for the ground state), 
or they constitute nonalternating divergent series with 
a subleading divergent alternating component (e.g., for 
ni = 0, 712 = 1, m = 0), or the series are alternating 
with a subleading divergent nonalternating component 
(e.g., for n\ = 1, n<x = 0, m = 0). The perturbation 
series for the Stark effect do not strictly fulfill the as- 
sumptions of Eq. (|l|), and the successful resummation 
of these series might indicate that the method introduced 
here is in fact more generally applicable. The large-order 
asymptotics of the perturbative coefficients for the Stark 
effect are given in Eqs. (4,5) in p3| . In quantum field 
theory, the alternating and nonalternating components 
correspond to ultraviolet (UV) and IR renormalons. Us- 
ing the first 20 coefficients of the perturbation series for 
the energy and evaluating the first 20 transforms accord- 
ing to Eq. (||), estimates for the real part of the energy 
(Stark energy shift) and the imaginary part of the en- 
ergy (decay width of the state) may be obtained. The 
apparent convergence of the first 20 transforms for the 
real part of the energy extends to 6-8 significant figures, 
whereas the convergence of the imaginary part is much 
slower (2-3 significant figures). In all cases considered, 
both the real and the imaginary part of the energy ob- 
tained by resummation compare favorably with values 
for the decay width obtained by numerical diagonaliza- 
tion of the Hamiltonian matrix [p~9| — ^T]| - Here we con- 
centrate on the decay width, the full calculation will be 
described in detail elsewhere. The atomic unit system is 
used in the sequel, as is customary for this type of cal- 
culation p^-|2l]|. In the atomic unit system, the unit of 
energy is or m c c 2 = 27.211 eV where a is the fine struc- 
ture constant, and the unit for the electric field is the field 
strength felt by an electron at a distance of one Bohr ra- 
dius dBohr to a nucleus of elementary charge, which is 
1/(4 7re ) (e/o| ohr ) = 5.142 x 10 n V/m (here, e is the 
permittivity of the vacuum) . 

Evaluations have been performed for all atomic lev- 
els and field strengths of Table III in p3|. Three ex- 
amples are presented here. For the ground state, at 
an electric field strength of E = 0.1 in atomic units, 
the imaginary part of the first 20 transforms calcu- 
lated according to Eq. (H) exhibits apparent convergence 
to r = 1.46(5) x 10~ which has to be compared to 
r = 1.45 x 10~ 2 obtained from numerical diagonaliza- 
tion of the Hamiltonian matrix JlSfl . For the L shell 
state with quantum numbers n\ = 0, n.2 = 1) Tn = 0, 
at a field strength of E = 0.004, the first 20 trans- 
forms exhibit apparent convergence to an imaginary part 
of r = 4.46(5) x 10~ 6 which compares favorably to 
r = 4.45 x 10~ 6 from pi] . The most interesting case 
is the state n\ = 1, TI2 = 0, m = 0, for which the non- 
alternating component of the perturbation series is sub- 
leading. At E = 0.006, resummation of the complete per- 



4 



turbation series (including the leading alternating part) 
leads to a decay width of T = 6.08(5) x 10~ 5 , which 
is again consistent with the result of T = 6.09 x 10~ 5 
from Mj. The contour C+i is crucial, due to poles lying 
off the real axis. 

With the help of Carleman's Theorem [Q it is possi- 
ble to formulate a criterion which guarantees that there 
is a one-to-one correspondence between a function and 
its associated asymptotic series (see for example p5| , 
Theorems XII. 17 and XII. 18 and the definition on p. 43 
in p(|, p. 410 in [^7j, or the comprehensive and eluci- 
dating review [[|). Let f(z) be a function which is ana- 
lytic in the interior and continuous on a sectorial region 
<S = {z\ \ arg(z)| < 7r/2 + e, < \z\ < R} of the complex 
plane for some e > 0. Let the function / have an asymp- 
totic expansion f(z) ~ S^Lo c ™ z " (f° r z ~ * The 
function / obeys a strong asymptotic condition (of order 
k) if there are suitable positive constants C and a such 
that \f{z) - £™ =0 c„z"| < Ca m + 1 [k(m+l)]\\z\ m + 1 
holds for all m and for all z E S. The validity of such a 
condition implies that the function f(z) is uniquely de- 
termined by its asymptotic series (see Theorem XII. 19 
of p6||). Typically, series which entail nonperturbative 
(imaginary) contributions do not fulfill the Carleman 
condition. The resulting ambiguity is reflected in the 
three integration contours in Fig. only one of which 
gives the physically correct result. 

It has not escaped our attention that specialized vari- 
ants of the method introduced here can be constructed in 
those cases where additional information about the per- 
turbative coefficients (large-order asymptotics, location 
of poles in the Borel plane, etc.) is available. 

Finite and consistent answers in quantum field theory 
are obtained after regularization, renormalization and re- 
summation. Using a resummation method, as shown in 
the Tables | and |lj, it is possible to go beyond the ac- 
curacy obtainable by optimal truncation of the pertur- 
bation series. The purpose of resummation is to even- 
tually reconstruct the full nonperturbative result from 
the divergent perturbation series (see also ||]). I have 
examined two physical examples, the QED effective ac- 
tion in a constant background electric field [Eq. (||)] and 
the Stark energy shift. The perturbation series for the 
Stark effect contains nonalternating and alternating di- 
vergent contributions, which correspond in their mathe- 
matical structure to IR and UV renormalons in quan- 
tum field theory, respectively. It has been shown in 
each case that complete nonperturbative results, includ- 
ing the pair-production amplitude for electron-positron 
pairs and the atomic decay width, can be inferred from 
the divergent nonalternating perturbation series by the 
resummation method defined in Eqs. (]]])-(§])• A mathe- 
matical model series (10), which simulates the expected 
large-order growth of perturbative coefficients in quan- 
tum field theory [see Eq. (Q)], can also be resummed by 
the proposed method (see Table |J) • In all cases consid- 
ered, the full nonperturbative result involves an imagi- 



nary part, whereas the perturbative coefficients are real. 
The advocated method of resummation makes use of the 
Pade approximation applied to the Borel transform of 
the divergent perturbation series. Advantage is taken of 
the special integration contours Cj (j = —1,0,1) shown 
in Fig. 0. The author acknowledges helpful discussions 
with G. Soff, P. J. Mohr and E. J. Weniger. 
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